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We show that pairs of widely separated interferometers are advantageous for measuring the 
Stokes parameter V oi & stochastic background of gravitational waves. This parameter charac- 
terizes asymmetry of amplitudes of right- and left-handed waves and generation of the asymmetry 
is closely related to parity violation in the early universe. The advantageous pairs include LIGO 
(Livingston)-LCGT and AIGO- Virgo that are relatively insensitive to flew (the simple intensity of 

' the background). Using at least three detectors, information of the intensity f^cw and the degree 

' of asymmetry V can be separately measured. 



I. INTRODUCTION 

Stochastic background of gravitational waves is one of the most important targets for gravitational wave astron- 
omy. In the last decade, the detection threshold for the background has been rapidly improved around ~ lOOHz by 
continuous upgrades of ground-based interferometers [l| . This trend will be continued with advent of next-generation 
0^ interferometers currently planned worldwide, such as advanced LIGO 2] and LCGT Due to the weakness of 
Q gravitational interaction, our universe is transparent to the background up to very early epoch, and we might uncover 
5^ interesting nature of the universe at extremely high-energy scales, through observational studies of the stochastic 
background. To extract the information as much as possible, we need to characterize the background efficiently in a 
. model independent manner, and investigation beyond simple spectral analysis might yield a great discovery. In this 
respect, circular polarization degree, which describes the asymmetry between the amplitudes of right- and left-handed 
I \ waves, may be a fundamental characteristic of the background to probe the early universe. Because the parity trans- 
^ • formation relates these two polarization modes, the asymmetry in the stochastic gravitational waves directly reflects 
\ a parity violation in the early universe, for instance, generated through the gravitational Chern-Simons term (e.g., 
CO ■ Q). In other words, one can detect a signature of parity violation by measuring the circular polarization degree 
m ■ of a gravitational wave background. Since the observed universe is highly isotropic and homogeneous, we shall fo- 
, cus on the monopole component of the circular polarization as our primary target, and report principle aspects for 
t its measurement with a network of ground-based interferometers (see [5], |6| for CMB polarization and [7] for space 

missions). 

o 



II. CIRCULAR POLARIZATION 

Let us first describe circular polarization of a gravitational wave background. We use a plane wave expansion of 



d \ the background as [8|, |9| 

h,,{t,x)= p df j^^dnhp{f,n)e^-^f'^-'+'^'^'^eP{n). (1) 



P=+,x ■ 



Here, the amplitude hp \s the mode coefficient that is stochastic and random variable. The bases for transverse- 
traceless tensor (P = -f , x ) are given as e+ — eg eg — and = eg -\- e^ eg with unit vectors eg 

and 60. These vectors are normal to the propagation direction n, associated with a right-handed Cartesian coordinate 
as usual. As an alternative characterization, we can use the circular polarization bases = (e"*" + ie^)/\/2 (right- 
handed mode) and = (e+ — ie^)/^/2 (left-handed mode) for the plane wave expansion ([1]). The corresponding 
amplitudes hp and hp are given by hp — {h+ — i /ix)/\/2 and hp — {h+ -f i hx)/V2- The ensemble average of their 
amplitudes is classified as 



{hp{f, n)hp{f, n'Y) \ _ Sn.n'Sfj' ( I{f, n) + V{f, n) 
{hp{f,n)hL{r,n'Y) ~ 47r I /(/, n) - T/(/, n) 



(2) 
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with the functions 6y,z being delta functions. In the above expression, the real function V characterizes the asym- 
metry between the amplitudes of right- and the left-handed waves, while the function /(> \V\) represents their total 
amplitude. Note that the other combinations such as {hjih*j^) and {hLh*f^) describe the linear polarization mode and 
are proportional to Q ± i [/, which constitute the well-known Stokes parameter, together with the /- and V- modes 
(see e.g., pX)| for electromagnetic counterpart). In this paper, we do not study the linear polarization Q ±iU , since 
they do not have an isotropic component. We will focus on the detectability of the isotropic components /(/) and 
V{f) as our primary target. Using the normalized logarithmic energy density of the background ^cwif) the 
two functions / and V are expressed as 

m = f^ow(/), vif) = i)ow(/) n(/), (3) 

where pc is the critical density of the Universe, pc — SHq/Stt with Hq = 70ft-7o km/sec/Mpc being the Hubble 
parameter. The ratio n(/) = V{f)/I{f) characterizes the circular polarization degree. For simplicity, we assume 
the flat spectra, flcwif) oc f'^ and n(/) cx f'^ as our fiducial model. Thus, our main interest is the simultaneous 
determinations or constraints on the parameters f^cw and 11. 

We next consider how to detect the isotropic components of /- and y-modes with laser interferome- 
ters. Let us recall that the output signal Sa of a detector a at the position Xa is written as Sa(/) — 
J2p=+ X Is^ ^^pif' ('^j /) 6* ^ir/n ajo ^ jjcre, the function is the beam pattern function and it represents 

the response of the detector to a polarization mode . Provided the data streams Sa and Sh taken from two detectors 
a and b, the detection of stochastic signals can be achieved by taking a cross-correlation, (sa(/)s&(/')*) = Cab{f)5fj' ■ 
Keeping the signals from the isotropic components, the correlation signal Cabif) is written as 

Cabif) = llMfW) + lV,abif)Vif), (4) 

where the quantity 7/ is the overlap function given by [1, Q 

li.abif) = ^ f dn [{F+Fr + F^Fr} e^^^^] , (5) 
»7r Js2 

with y = 2-Kf D/c. Here, we have expressed Dm {D: distance, m: unit vector). Similarly, the function 

lv,ab{f) is obtained by replacing the kernel [• • •] in Eq. © with [i {F+F^* - F^ F+*} e'v'^'^] . 



III. OVERLAP FUNCTIONS FOR GROUND BASED DETECTORS 



Now, specifically consider the response of an L-shaped interferometer a on the Earth. We assume that the detector 
has two orthogonal arms with equal arm- length. Denoting the unit vectors parallel to the two arms by u and v, 
the beam pattern function takes a simple form as F^ = da : e^(n) with da = {u (S) u — v ^ ■i')/2, where the colon 
represents a double contraction. This expression is always valid as long as the wavelength of the gravitational waves 
for our interest is much longer than the arm-length of the detectors. In this paper we study the following five ongoing 
(and planned) kilometer-size interferometers as concrete examples; AIGO(A), LCGT(C), LIGO-Hanford(H), LIGO- 
Livingston(L) and Virgo(V) (see e.g. [Tlj for their basic information). Hereafter, we mainly use their abbreviations 
(A,C,H,L,V). 

For the isotropic component of the stochastic background, only the relative configuration of two detectors is relevant 
with the correlation signal Cab and we do not care about the overall rotation. Hence, the sensitivity of each pair 
of detectors to the stochastic background can be characterized by the three angular parameters {f3,ai,a2) shown in 
Fig. [1] Here, /3 is the separation angle between two detectors measured from the center of the Earth. The angle cti 
(tT2) is the orientation of the bisector of two arms for detector a (6) measured in counter-clockwise manner relative 
to the great circle connecting a and b. Their distance is given by D — 2Re sm{f3/2) {Re — 6400km : the radius of 
the Earth), which determines a characteristic frequency fo = 0/(271 D) for the overlap functions. Following Ref.jl], 
we define the angles 

A EE (fxi + ^2)/2, 5={ai- a2)/2. (6) 

The geometrical information about pairs of detectors among the five interferometers is presented in Table IH 

In the expression ([5]), the angular integral can be performed analytically with explicit forms of the pattern functions. 
A long but straightforward calculation leads to Q 



li,ab = Qi{y,f3) cos(4,5) + e2(y,/3) cos(4A), 



(7) 



FIG. 1: Geometrical configuration of ground-based detectors a and b for the cross-correlation analysis. Detector planes are 
tangential to the Earth. Two detectors a and 6 are separated by the angle j3 measured from the center of the Earth. The 
angles cti and (T2 describe the orientation of bisectors of interferometers in a counter-clockwise manner relative to the great 
circle joining two sites. 
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TABLE 1: Upper right cos(4(5)). Lower left (cos(4A), sin(4A)). 



with ei(y,/3) = cos* (I) (jo + fj2 + 112^4) , and Q2iv.fi) = (-fjo + %k " M^) + il^ - 1^ - Mi^) ^os/? 
-|- — — gfg ^4) cos(2/3). The function j„ is the n-th spherical Bessel function with its argument y = f / fn- 

[3] • On the other hand, the overlap function for the F-mode is given by 

7v,a6 = e3(2/,/3) sin(4A) (8) 

with 63(7/, (3) = — sin \_{^h + fja) + (ii + f Ja) cos 0\ . In Fig. [21 the overlap functions for the two representative 

pairs are shown in top (HL) and middle (CL) panels. 

Here, we give a simple interpretation for the angular dependence of Eqs. ^ and ([5]). The beam pattern functions 
and are given by linear combinations of ( cos(2(Ti), sin(2(Ti) ) and ( cos(2cr2), sin(2(T2) ) respectively, reflecting 
their spin-2 like nature. Then, with Eq. ([5]) and addition formulas of trigonometric functions, the overlap functions 
should be linear combinations of cos[2(cri ± (T2)] and sin[2((Ti ± (T2)], namely, cos(4A), cos(4(5), sin(4A) and sin(4(5). 
Since the expectation value Cabif) is a real function for our beam pattern functions, we have (sasj) — (sbsl). This 
essentially results in replacing the roles of ai and (T2 , and the functions 7/ and cannot contain terms proportional 
to sin(4^) — sin[2(cri — 0-2)]. 

On the other hand, while the observable Cabif) and the amplitude / are invariant under the parity transformation of 
a coordinate system, the sign of the parameter V flips, because the transformation interchanges right-and left-handed 
waves. Therefore, the function ^v.ab must change its sign while keeping the quantity Cabif) invariant. Geometrically, 
this corresponds to the re-definition of the azimuthal angles cri 2 in a clockwise direction (or putting ai — *■ — Ui and 
(72 — 172). As a result, the function "fv,ab should be odd functions of S and A, and it must be proportional to 
sin(4A) as shown in Eq. ([S]) (the term proportional to sin(4(5) is already prohibited as explained earlier). With similar 
arguments, we find that the function 7/ is a linear combination of cos(4A) and cos(4(5) as in Eq. ([7]). 

A. Special cases 

To stress the importance of the geometric configuration, it is instructive to consider several simple examples for 
idealistic pair of detectors. When a pair of detectors is co-located (/3 = 0° and D = 0), the functions (81, Q2) defined 
after Eq. ([7]) become (1,0) and we have ji^ab — cos(4(5). The identity 82 = at /3 = 0° imphes that the function 
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FIG. 2: Overlap functions for the un-polarized I mode (dashed curves), and the circularly polarized V-mode (solid curves). The 
upper panel shows the results for the Hanford-Livingston (HL) pair (the characteristic frequency fo ~ lOOHz). The middle one 
is results for the LCGT-Livingston (CL) pair (/d — 31Hz). The normalized SNRs Si^v (with the adv LIGO noise spectrum) 
are also presented. The bottom one show the compiled functions F/.v (ea. (|10|l ') made from both pairs. 

7/ depends very weakly on the parameter A at a small angle /3 ~ 0°. On the other hand, the overlap function jv.ab 
always vanishes for a pair of detectors in the same plane {(3 = 0). This is even true with a finite separation D ^ 0. 
This exact cancellation comes from the geometric symmetry of the beam pattern function with respect to the detector 
plane 0,111. 

For two detectors at antipodal positions (/3 = 180°), we have 6i = and the angle S becomes geometrically 
meaningless. One can expect that the function ji^ab is almost proportional to cos(4A) near f3 = 180°. 

B. Broadband SNR 

Now, we turn to focus on a broadband sensitivity to the /- and T^-modes. In the weak signal limit, the total 
signal-to-noise ratio (SNR) for the correlation signal Cab{f) is given by Q 

with X = 7/ r^Gw + 7y f^GwII. Wc denote the noise spectra for detectors a and b by Na{f) and Ni,{f), assuming 
no noise correlation between them. In what follows, for simplicity of our analysis, we further assume that all the 
detectors have the same sensitivity comparable to the noise spectral curves of advanced LIGO. The analytical fit 
from Fig. 1 of Ref.[i] leads to N{f) lO^*"^ (//lOHz)"* + lO^^'^ ^s (//lO^Hz)^^'^ Hz^i for lOHz < / < 240 Hz, 

N{f) = lO"'*^ (//lO^Hz)^ Hz-i for 240 Hz < / < 3, 000 Hz, and otherwise N{f) = oo. Note that the combination 
f^N{f)^ becomes minimum around / ~ 50Hz with its bandwidth A/ ^ lOOHz. For a pair of coincident detectors 
(i.e., ji^ab = 1 and jv,ab = 0), the total SNR is evaluated by setting X = ficw in Eq. ([H]), and we numerically obtain 
SNRo = 4.8 (Tobs/Syr)'/' (n^^^JlQ-^) . 

The total SNR depends strongly on model parameters of the background, including the polarization degree H. In 
order to present our numerical results concisely, we first calculate SNRjj ^jf, by plugging X = J{iy}^ab into Eq. ^ 
and then normalize them as S[j,v}.ab = SNR|/ Qf,/SNRo. The normalized SNRs can be regarded as rms values of 
l{i,v}.ab with a weight function [f^ N{f )'^]^^ . 

C. Optimal configuration 

Let us discuss optimal configurations of two detectors (a, 6) for measuring the /- and y-modes with the correlation 
signal Cab- There are two relevant issues here: maximization of the signals Si^ab and Sv.ab, and switching off either of 
them {Si^ab = or Si^ab = 0) for their decomposition. To deal with the situation comprehensively, we consider how to 
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set the second detector b relative to the fixed first one a with a given separation angle p. In this case, the sensitivity 
to the /- and V-modes is characterized by the remaining adjustable parameters, CTi and (72 • The former determines 
the position of the detector b, while the latter specifies its orientation (see Fig. [1]). Based on the expressions ^ and 
([5]), one finds that there are three possibilities for the optimal detector orientation: cos(4A) = — cos(4(5) — ±1 (type 
I) or cos(4A) = cos(4(5) = ±1 (type II) to maximize the normalized SNR Sj^ab [1], and cos (4A) = cos (4(5) = (type 
III) to erase the contribution fi-om /-mode. For type I, the solutions of the two angles ai.2 are cti = (72 = 45° (mod 
90°) and the detector b must be sited in one of the two great circles passing through the detector a, parallel to one 
of the two arms. For type II, the second detector must reside in two great circles parallel or perpendicular to the 
bisecting line of each detector. Similarly, the type III configuration is realized by placing the second detector on one 
of the four great circles defined for types I and II, with rotating 45° relative to the first detector. 

Note that the sensitivity to the T^-mode is automatically switched off for the type I and II configurations and is 
conversely maximized for the type III configuration. This is because the normalized SNR Sv,ab is proportional to 
sin(4A). While a definite detection of a weak V-mode signal requires a careful removal of the /-mode signal from 
observed data, it turns out that the geometrical requirement for type III configuration is severe. As we see later, 
however, we can easily control the contribution from the /- (or V-)mode by introducing a third detector. 

In Fig. [31 we present the normalized SNRs for the optimal geometries; types I, II and III (short-dashed, long-dashed, 
and solid curves, respectively). One noticeable point is that a widely separated (/? ~ 180°) pair is powerful to search 
for the V-mode (recall the cancellation 7v = at /3 = 0). To reduce the contribution from the /-mode, pairs that are 
usually disadvantageous to measuring the total intensity J^gw now play a very important role. In Fig. [31 we also show 
the normalized SNRs for representative pairs made from the five detectors, in which several interesting combinations 
are found. The HL (with cos(4(5) ~ 1 and sin(4A) ~ 0.93) realizes nearly maximum values simultaneously for Sj^ab 
and Sv,ab at its separation f3 = 27.2°. This is because Si^at is mainly determined by the angle i5 at a small /3, while 
Sv,ab depends only on A. The CL has good sensitivity to the T^-mode and relatively insensitive to the /-mode with 
sin(4A) ^ 1. In contrast, AH is almost insensitive to the V mode with sin4A = —0.007. In this sense, LCGT and 
AIGO detectors are suitably oriented to probe the /- and V-modes, respectively. 




FIG. 3: Normalized signal to noise ratios {Si^ab and Sv,ab) with optimal configurations for the 7-mode (short dashed curve: 
type I, long dashed curve: type II) and for the V-mode (solid curve: type III with setting IT = 1 for illustrative purpose). We 
use the noise curve for the advanced LIGO. For each detector pair, Si and 5*^ are given with a triangle and a circle respectively 
at its separation /?. There are four other pairs not shown here; CH with [Si, Sv) = (0.04, 0.08), LV with (0.08,0.04), HV with 
(0.07,0.06) and CV with (0.09,0.04). 



D. Separating /- and 1^-modes 



As a final mention, we will address the issue of /- and y-mode separation by combining several pairs of detectors. 
For preliminary investigation, we consider the case that two pairs of interferometers (a, b) and (c, d) are available. 
Detectors a and c can be identical, but we need at least three independent detectors for the study below. First note that 
the correlation signals are given by Coh(/) = li.abif) I{f)+lv,ab{f) V{f) and Ccd{f) = li,cd{f) I{f)+lv,cd{f) ^(/)- 
From this, one can easily find that the contribution from the /-mode is canceled by taking a combination W = 
"yi,abCcd — Ji,cdCab = {Tv,cdji,ab — Jv,abli.cd) V{f). Thc Statistical analysis based on the combination W would be 
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a robust approach for actual V-mode search, although a further refinement may be possible by combining more pairs, 
which we will report elsewhere. 

Since the rms amplitude of the detector noise for the combination W becomes N{f) (7I + cdV^^^ define the 
compiled overlap function for the V mode by 

i V,ab:cd = r 2 , ^2 11/2 ' (10) 

This expression should be used in Eq. ([9|) when evaluating the broadband SNR for the V-mode with the combination 
W . In a similar way, we define the compiled function T i^ab:cd for the / mode by interchanging the subscripts V 
and / in Eq. (|10p . Bottom panel of Fig. [2] shows the compiled overlap functions ^ {iy}.ab:cd from two pairs of 
detectors, CL-HL. With this combination, the normalized SNR becomes 0.11 for the V-mode and 0.31 for the I- 
mode. Using numerical results below eq.([5]), the detection limit for the polarization degree 11 is given as 11 = 
(T/3yr)~^/^(S'A^_Ry/5)(riGw/i7o/10^^)^^ with signal-to-noise ratio SNRy- These numerical results are almost the 
same values as in Sy for CL and 5*/ for HL, and in this sense, the /-, V-mode separation can be performed efficiently 
with naively expected sensitivities S^jy},ab- Note that the other combinations, such as AV-HL, AV-HV and CL-HV, 
also provide the normalized value ^ 0.11 for the V-mode, but AH-AL has only 0.015. 

In summary, we reported principle aspects for measuring a circular polarization degree of a gravitational wave 
background that is related to parity violation. We find that pairs of ground-based interferometers that are widely 
separated and relatively insensitive to the total intensity f^cw are advantageous for the measurement. With at least 
three detectors, the polarization degree and the intensity ficw can be separately detected. 
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